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iy-j , We obtain correction terms to the large N asymptotic expansions of the eigenvalue density for 

■ the Gaussian unitary and Laguerre unitary ensembles of random N x N matrices, both in the bulk 

of the spectrum and near the spectral edge. This is achieved by using the well known orthogonal 
£SJ , polynomial expression for the kernel to construct a double contour integral representation for the 

density, to which we apply the saddle point method. The main correction to the bulk density is 
oscillatory in N and depends on the distribution function of the limiting density, while the corrections 
^r* to the Airy kernel at the soft edge are again expressed in terms of the Airy function and its first 

derivative. We demonstrate numerically that these expansions are very accurate. A matching is 
\& • exhibited between the asymptotic expansion of the bulk density, expanded about the edge, and the 

asymptotic expansion of the edge density, expanded into the bulk. 

PACS numbers: 02.50.Cw,05.90.+m,02.30.Gp 
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We consider in this paper two classical ensembles of random matrices, the Gaussian unitary ensemble (GUE), and 
the Laguerre unitary ensemble (LUE). These ensembles can be characterized by their joint eigenvalue probability 
' density functions 
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P N (xi,...,x N ) <xY\_un(xi) I j (x k -Xj) 2 , xiGtt, (1) 

1=1 l<j<k<N 



Jexp(-2iVx 2 ), GUE, 
ujn(x) = < (2) 
v ' \x a exp(-4A^.T), LUE, v ' 



n= R ' ^ GUE ' (3) 

I (0,oo), LUE. v ; 

The GUE consists of N x N Hermitian matrices with independent normally distributed entries on and above the 
diagonal. It is the cornerstone of random matrix theorjiSi^. The LUE has fundamental applications in mathematical 
statistics and quantum field theory since it includes Wishart matrices and the Chiral GUE as special cases (the latter 
after a straightforward change of variables); see e.g. Refi 

We are interested in the large iV behavior of the marginal eigenvalue probability density pm{x), which we hereafter 
refer to simply as "the density" , and which is defined by 



p N {x) := / P N (x,X2,---,x N )dx2.--dx N . (4) 



The function N pn{x) can be interpreted as the number density of eigenvalues near the point x. We also remark 
that for the GUE, TV pn(x) is equal to the number density of a harmonically trapped system of either non- interacting 
fermions or impenetrable bosons^. There is a similar interpretation for the LUE in terms of a Calogero-Sutherland 
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type model^. For recent advances in asymptotic questions related to these interpretations, complementary to the 
present study, see Refsi 6 i 7 i 8 

As background to the present study we note that aspects of the large N form of pn(x) first arose in studies of field 
theories related to Hermitian matrix models^. There, for the GUE the large N asymptotic expansion of the moments 

niN(p) ■= / x p p N (x)dx (p=l,2, ...) 
Jn 

was sort. By a graphical expansion of the matrix integral, involving cataloging the corresponding maps according to 
their genus, it was predicted that for certain coefficients a,2j{p), 

p/2 , , 

< UE b) = E^T> (P = 2,4,...) (5) 

3=0 

(the odd moments of course vanish). Analogous considerations in the case of the LUE 10 show that 

^ UE ( P ) = E%^ (6) 

for certain coefficients dj{p,a). Observe in particular that (JHJ) contains only even inverse powers in N, while © 
contains both even and odd inverse powers in N. 

The graphical methods allow ao(p) in (JSJ and ao(p,a) in JSJ to be computed in terms of binomial coefficients for 
all p = 0, 1, . . . . This knowledge in turn can used (see e.gi) to prove that in the limit N — ► oo and with x fixed 



1,1], GUE, 




^) : =J 1 T o Pjv(:E) = i-A/--l ,3 6(0,1], LUE, ( ? ) 

| 7T V X 

otherwise. 

The first functional form in (JJJ is referred to as the Wigner semi-circle law, while the second is sometimes named 
after Marcenko-Pastur. See e.g. RefsiS* 3 .. 

The expansions (J5J and © provide a motivation to undertake a study of the asymptotic form of (@J . In the case 
of the GUE such a result has been given by Kalish and Braak 11 . It states that for |x| < 1 and fixed 

, s 2cosf2A^7rF(a;)l 1 / 1 \ 

where 

x 1 
P(x) = 1 + — p(x) Arccos(a;). 

2 7T 

Thus one sees that unlike the situation with the moments (0, the leading correction term is 0(1/N). Of course this 
term is oscillatory so one might anticipate that after integration it contributes at a higher order. However inspection 
of (JSJ) reveals that the situation is more complex: the oscillatory term is not integrable at the endpoints of the support 
|z| = 1. Indeed, it is well known that with the boundary of the eigenvalue support taken as the origin, a scaling 
regime distinct from that of the bulk becomes relevant. Explicitly, with Ai(x) denoting the Airy function, it has been 
proved thalii^ 

Nl/3 „gue L e \ _ l; „ (2ao 1/3 .UE {, , e 



lim Pn 1 H ^7 = lim Pnx- , , > 

n^oo 2 HN V 2Ar 2 /3 J w^oo 2 HN V (2A) 2 /3 J (9) 

= [Ai'(e)] 2 -e[Ai(C)] 2 , 

where ^ is fixed. In view of the breakdown of JHj in the vicinity of the spectrum edge, (referred to as the soft edge, 
since although it defines the edge of the support of p(x), for any finite N there is a nonzero probability of finding 
eigenvalues lying beyond it), we are thus led to also investigate the large N asymptotic expansion extending the limit 
law ©. 
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At a technical level, the main achievement of this paper is the derivation of the first correction terms to the limit 
laws (J7J and ©. We do this by utilizing the well known orthogonal polynomial expression for pn(x) to obtain a 
double integral representation which is amenable to the saddle point method. In the bulk, i.e. in the interior of the 
support of p(x), we show that the asymptotic series progresses in powers of 1/-/V, and we obtain the explicit form 
of the 1/A correction. We find for the LUE that the coefficient of the 1/A term consists of a component which 
is oscillatory in A as well a component which is non-oscillatory in A, whereas for the GUE it consists of only an 
oscillatory component, as shown in |J5J. For the soft edge we will see that the asymptotic series progresses in powers 
of A" 1 / 3 , and we obtain explicit expressions for the coefficients of the A and TV 2 / 3 terms, which again involve 
Airy functions. 

Due to the similarity in the structure of Pn(x) for the GUE and LUE, it is convenient to consider both cases 
simultaneously to avoid unnecessary repetition, and so at each step of our presentation we discuss the GUE and LUE 
in parallel. In Section^ we discuss the double contour integral expression for pn(x) to which we shall apply the 
saddle point method. Section llTTl contains our discussion of the asymptotics of pn(x) in the bulk while Section llVl 
discusses the soft edge. In Section we discuss the extent to which our expansions in the bulk match up with those 
for the soft edge. 



II. CONTOUR INTEGRAL EXPRESSION FOR p N (x) 

For the unitary ensembles there is a well known and very neat expression for pn{x) in terms of orthogonal poly- 
nomials, valid for any A and x € SI. If we let {^j(x)}j°^ denote the monic polynomials orthogonal with respect to 
u>n{x) on Q,, then 

P N ^ X ) = A n N ^ X \\2 ^ 7r N( x ) 7r N-i{x) - n' N _ l (x)ir N (x)]. (10) 

iV ||7TjV-l|r 

The norm in the denominator of l|10|l is just the L 2 -norm associated with wjv(x) and fl. For a derivation of i|10|) the 
reader is referred to Refsi2*£. We remark at this point that there is no universally agreed scale by which the GUE 
and LUE are defined. To match our choice of scale and notation in @ to that employed in Refi for instance, we 
observe that 



Pn(x) 



' V2N- 1 / 2 P N (V2Nx, VMx), GUE, 
AP N (4Nx,ANx), LUE, 



where Pn(x,u) (not to be confused with our definition Q above) is the kernel defined in Chapter 4 of RcfA (the 
kernel is often also denoted Kn(x,u) in the literature) . 

To investigate the large A behavior of pn{x) it is obviously advantageous to start with the expression ifTT))) rather 
than with the (A — l)-fold integral Q). The iTN+j-i(x) can be expressed in terms of the standard Hermite and 
Laguerre polynomials found in Szego's classic bookA^ as follows 

( 2 - 3 ( N +j-V/ 2 N-( N +j- 1 )/ 2 H N+j _ 1 (V2Nx), GUE, 
tjv+j-iW - y^N+j-i^ + j _ i)l(4j V )-JV-J+i i W i _ 1 (4JVx), LUE. (U) 

The required asymptotic expansions of the scaled Hermite and Laguerre polynomials appearing in l|llfl are known 
to any order both in the bulk and near the soft edg o 14 ' 15 , and such asymptotic expansions of scaled orthogonal 
polynomials are now generically said to be of Plancherel-Rotach type (Plancherel and Rotach were the first to compute 
such asymptotics for the Hermite polynomials). It is reasonable to assume that the most straightforward procedure 
to obtain the desired asymptotic corrections for pn{x) in each region of interest is to simply insert the corresponding 
asymptotic expansion for 7r/v+j-i(x) into 1|10[) . While this is certainly legitimate in principle, and does indeed recover 
the leading term fairly easily, to derive the correction terms it turns out that such a procedure is rather tedious, and 
provides little if any insight into the resulting expressions. The source of the complication is that the asymptotic 
expansions for contain a large amount of superfluous information which is canceled when the expansions 

are substituted into l|10fl . To avoid this, we shall pursue a related, but more direct route. 

The Plancherel-Rotach asymptotics for the Hermite and Laguerre polynomials were originally derived by first 
expressing the polynomials in terms of contour integrals, and then applying the saddle point method. By suitably 
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massaging the standard results in Szego's booisi^ one finds that 

(r tJz p~Nz 2 /2 
c ^ N) tZl e2N -^r^ GUE, 
r dz (z + 2) N + a n n 3 '- 1 (12) 
I t\N+j-l„ <AT\ 1 aZ -2NZX V- Z > I 1 I M TTTP 

27Ti Z^+i 

(JV + j -1)! 
(2iV) 7 



'•-■■ v ' := v.;j^--(' - 



In both cases the contour of integration is a closed positively oriented contour which encircles the origin; in the 
Laguerre case we further demand that it not contain the point z — — 2. 

Instead of applying the saddle point method to l|12l) and then substituting the expansions into (|10|) , wc shall first 
insert the contour integrals (|12[) into (|10|l to obtain a double integral expression for p^(x), and then perform the 
saddle point method on this double integral. To highlight the similarity between the GUE and LUE it is convenient 
in the GUE case to substitute the contour integral for 7Tjv+j-i(— x) into lfTU|) rather than that for 7rjv + j_i (x); since 
Hj(—x) — (~iy Hj(x) this ruse is perfectly harmless. This results in 

p N {x) = 2 C °[ N)Cl \^ ] ] - u N {x)J N {x), (14) 

where 

Jn{x) := / 2Vi f ^e NSM+NSM G(^,z 2 ) (15) 

and 

f-2«-log(*)-zV2, GUE, 

K ' h \-2zx-log(z)+log(l + z/2), LUE, U ° J 

G(z 1 ,z 2 ):=u(z 1 )u(z 2 ) , (17) 

u(z):=l)> , GUE ' (18) 

V ; \(l + z/2) a -\ LUE. V ; 

The remainder of this paper will involve a careful asymptotic analysis of the double integral i|15|) ■ 

Before proceeding we note that it is straightforward to show, using standard results in the orthogonal polynomial 
literature^, that 

(o2AT-3/2 atN+1/2 
7= GUE > 
f yfn2N +a -i (19) 



] (4iV) 2 



t T(N) T(N + a) ' 
and hence the asymptotics of the prefactors in (114ft is 



LUE, 



} Cp(N} CyjN) 
' hN-l\\ 2 



-Ni~ N 



T(N), 
A N+a N a T(l + N) 



T(N + a) 



GUE, 



, LUE, 



2Ne 



-N 



4 JV+a N 



1 + 
1 - 



1 



o 



12 N ' ~ \N 2 
(a — 1) a 



1 



2N 



O 



N 2 



GUE, 
LUE. 



(20) 



(21) 
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Saddle points 

Before applying the saddle point method to l|15l) we need to identify and classify the saddle points of l|16fl . The 
functions S(z,x) in general have two saddle points at z = z± where 

f-zit !/(*), GUE, 
±- )-l±iu(x), LUE, 1 ' 

and 



VI -x 2 , GUE, 

"'" := < ' ■. LUE. <23) 



We note that for both the GUE and LUE we have 



v{x) = ^p(x), for |a;| < 1, (24) 

with p(x) as defined in Q. 
Since for \x\ < 1 we have 

■S"(z±,aO _ jv(x)e ±i ^- AlCBiD ^, GUE, 

2 ~ \2x 2 V {x)e ± ™' 2 , LUE, ( j 

the saddle points z — z± are both simple when |x| < 1, i.e. S"(z±,x) ^ 0. However, when x — 1 the two simple 
saddle points given in l|22|) coalesce to z± = — 1 and S"(z±, 1) vanishes, so we obtain one double saddle point in this 
case. Thus we already see why the regions |x| < 1 and x ~ 1 have qualitatively distinct asymptotic behavior. Simple 
saddle points generically produce Gaussian integrals whereas double saddle points generically produce Airy functions 
(see e.gji&). 

III. BULK ASYMPTOTICS FOR THE GUE AND LUE 

In the bulk of the spectrum, i.e. for \x\ < 1, we hold x fixed and investigate the asymptotics of (|15(l as N becomes 
large. From l|22|l we see that there are two distinct simple saddle points of S(z,x), which form a complex conjugate 
pair in this case. Let's define S± := S(z±,x). Then since 

Re(S+) = Re(SL), (26) 

both saddle points contribute to the same order and we deform our contour through both of them. Denoting by f2± 
a contour passing through z± along a path of steepest descent, the standard arguments of the saddle point method 
(see for e.g. Refii&) yield 



J N (X)=[I P-.+ ! p:)l[ p:+f ^ ) e * S M+ N S( Z2<X ) ^ ^ + 0{e 2NKe { S + )-N % (2?) 

y f dzi f ^ e NSiz 1 ^)+N8iz a ^) G ^ Z2)+0 ^lNB^S + )-Ne^ (2g) 

.rr 1 „js!„, 2lT Jn„„ 2 n 



for suitably small e > 0. 

We now need to parameterize the contours f2± . From (|25|) we see that if we set 

_ fArcsin(x)/2, GUE, 
~~ ] tt/4, LUE, 



(29) 



then a suitable parameterization of £l± is 

z = z ± + e ±ie t, t€[-r),r)}, (30) 
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for sufficiently small i] > 0. In (|29(l the function Arcsin : [—1, 1] — > [— tt/2, it/2] denotes the principle branch of arcsine. 
With the parameterization l|3U|) the contour f2 + is traversed in the negative direction, so we need to compensate for 
this with an explicit minus sign. 

By choosing rj sufficiently small S(z, x) is analytic on fl± and so using the parameterization l|3U[l we see that 

S(z,x) = S± - at 2 - at 2 tp±(t), z £ fi ±) (31) 

where 



v(x), GUE, 
2x 2 v{x) 1 LUE, 



(32) 



and 



-=2^sm^)M e 1 ■ (33) 

We note that (f-(t) — <p+(t), where * denotes complex conjugation. 
It is a straightforward exercise to show that S_ = S*L, and that 

Re(5+) = {2 +x2 ' GUE ' (34) 
[2x-log(2), LUE, 

lm(S+) = -nP(x), (35) 

where 

P(x) := [ p{t) dt, (36) 

Jx 

X 1 

l + -p(x) Arccos(x), GUE, 

1 + x p(x) Arccos(v / s), LUE. 

Here xq is the left edge of the support of p(x) given in J7J, i.e. xq = 0, 1, for the LUE and GUE respectively. We note 
that P{x) is the probability distribution function corresponding to p(x). The limiting distribution function P{x) will 
play a significant role in the bulk asymptotic expansion of pjv(x). 
With the results l|31|) and l|35|l for S(z,x), and the definitions 

G„{t) := G(zi,z 2 ) . 1 ^ 1+e i"ie tl , (38) 

/V -aJVt? 
— — dt u (39) 
-n 2vr 

we can make the change of variables l|30|l in (|28|l to obtain 

J N ( X ) = e ^R C (S + ) E (B) E (B) J2 - ai a 2 e ^i+-2)(.e-N,P( x )) e -aNtl Vai ( tl )-aNtl^ 2 (t 2 ) Gait) 

ae{+i,-i} 2 (40) 

A small amount of massaging shows that the — , + term in l|40|) is the complex conjugate of the +, — term, and likewise 
the — , — term is conjugate to the +, + term. Suppose now that we define the function F a (X,t) by 

F a (X,t) := e x ^"^e X2V '^G a {t) (41) 

then 

J N (X) =2e 2NR < S +~> ReE[ B) Ef } ( F + _(A,t) - e H0-N,P(x)) p+ +M \ + Q^iVRe^-eiVs 

[ ' \ i = -aNt 2 ' \i = -aNt 2 J 

(42) 
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We shall discuss the purpose of the parameters Ai shortly. Note that the two terms in (|42|l are qualitatively distinct 
- the second term is oscillatory in N while the first is not. 

It is convenient to pause for a moment and multiply 142JI by the the explicit forms for the prefactors required in 
(|14H using and l|34fl to obtain the corresponding expression for pn{x). Since 



lFJv-iir 



N 



N 2 



with 




we have 
p N {x) = h 



1 



N 

0(e^ N ). 



O 



1 

N 2 



hi = 



2ReiVE< t B) E 2 £i) <j F+,_(A,i) 



(43) 

(44) 
(45) 



2i(8-NirP(x)) 



(46) 



The introduction of the auxiliary variables Ai, A2 in 1|41[) is a common ruse applied in the saddle point method (see 
e.g. Ref^&) which we now discuss. Suppose that we construct the Maclaurin expansion in ti,t 2 of F a (X,t) with A 
considered as a fixed parameter 



k J-k 



j=0 k=0 w J 



Qj k 

R^P^ (A ' S1 ' S2) 



Si ,S2=0 



•oc*? 1 «?; 



Pl+P2=P+l 



(47) 



If we now set Ai = —aNtf in Ij47(l and perform the integrations required in l|46|l then we find that each term 
corresponding to a given value of j in l|47|) has the same resulting N dependence. This then gives a systematic way of 
obtaining the corrections out to any given order in N. To see why this occurs, first note that we need only consider 
the terms in Ij47(l for which both j and k are even since any odd monomials are annihilated by H39|l , and then further 
note that with A, = —aNt 2 we have 



NE\V h 2mi A, 1 An 2 = (-l) h+h F(mi +h + y 2 ^ m2 +h± W N -( mi+m2 ) Q( -y N) 

I z 1 z 1 z \ J ^ 7r 2 a mi+m 2 + l v " 



(48) 



for any € N. Hence despite the fact that various powers of Ai and A2 arise when F a (\,t) is differentiated, 

for a given value of j all terms end up with the same N dependence after setting Ai = —aNt 2 and performing the 
integrations. 

Hence, substituting (|47|) into l|4tj[l one can construct the asymptotic series for pn{x) out to any desired order. We 
shall explicitly construct this series out to order 1/N but the generalization to higher orders is obvious. However, as 
we shall see the resulting asymptotic expansions obtained by keeping only the 1/N correction are already extremely 
good, and it's quite likely that optimal truncation occurs at this order, providing little incentive to construct higher 
order corrections. 



Let us denote by Cm (x) the coefficient of 1/N m in the 1 /N expansion generated by acting on l|47|) with N E^"' E. 



Then since the only terms which contribute are those for which j and k are even we have 



Cm (3-) 

TV" 1 



1 



£j(2fc)!(2[m-fc])! 



NE^E^tf tf m - k) 



Qlk Ql(m — k) 



dS! 2k ds 2 2{ - m - k) 



(49) 



\ i= -aN tj 



We can re-express pn (x) from (|46|l in terms of the c~m (x) as follows 



PN (x) =2h Rc { C ^- 1} (ar)} - 2 h Re ^-n.p^)) c (M) (a;) j 



2h, 



oRe{ 



e 2i(8-NirP(x)) c , 



1 

N 



(50) 
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-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 



FIG. 1: Comparison of the asymptotic expansion l|f>ip. shown as the dashed line, and the exact result il l (II . shown as the solid 
line, for the eigenvalue density of the GUE with N = 10. 



It is not hard to show that the term 2 ho Re j^Cp 1 ' ^(x)l is equal to the limiting density p(x) from J2} when |a;| < 1 

as required, and also that Cq (x) vanishes identically (this is actually obvious from (|17|l and l|38j l). The explicit 
construction of the remaining ct' (x) required in (|50|l is straightforward and we finally obtain the following. 



Proposition 1. Let Pn{x) be as defined in ^j), and let x be fixed with \x\ < 1. Then as N 
For the GUE 



p N (x) = p(x) 



2cos[2 7Vtt P(x)] 1 
tt 3 p 2 (x) N 



0[ N* 



oo we have the following: 



(51) 



while for the LUE 

p N (x) = p(x) 



cos(2 N 7T P{x) - a tt [1 + x p(x) - P(x)]) 

7T 3 X 2 p 2 (x) 

cos(2 NttP(x) — 2 a Arccos(y / a;)) a 



tt 2 x p(x) J N 



1 

N 2 



1 

N 



O 



1 



(52) 



7T 3 x 2 p 2 (x) 4(1 — x) 

where p(x) is given in and P(x) is the corresponding probability distribution function, given explicitly in 

As remarked in the Introduction, the result 1)51(1 was obtained previously in 11 again by steepest descent, but starting 
from an integral representation derived by super-symmetric arguments rather than orthogonal polynomials. 

It is interesting to note that it is the distribution function P(x) of the limiting density p{x) which controls the 
large N oscillations in the 1/N correction to pn(x). Note also that while the non-oscillatory correction vanishes at 
order 1/N for the GUE leaving only an oscillatory correction at this order, the LUE has both an oscillatory and a 
non-oscillatory component to its 1 /N correction. 

To demonstrate to the reader just how accurate the expansions given by Proposition^ are we provide in Figures 
n and |21 a numerical comparison of the asymptotic expansions with the exact results computed using the expression 
HI in terms of orthogonal polynomials . 
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_0 5 I i i i i i i i i i I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x 

FIG. 2: Comparison of the asymptotic expansion (1")2> . shown as the dashed line, and the exact result il l (II . shown as the solid 
line, for the eigenvalue density of the LUE with a = 1/2 and N — 10. 

IV. SOFT EDGE ASYMPTOTICS FOR THE GUE AND LUE 

The appropriate scaling to elucidate the behavior of Pn(%) near the soft edge is to set x = 1 + £/iV 2 / 3 for fixed £, 
as appears in (jSJ). Substituting such a scaling into 1(16(1 we find 

NsfzA + ^Jfc) =-Cb 1/3 N^z + NS(z), (53) 

where we've introduced the shorthand S(z) := 5(2,1). Here b > is a free parameter that we can fix later as 
convenient. Defining 

Jn(0--=Jn{1 + ^j^), (54) 

we see that (153(1 leads to 

J N {i)= i^-.^v{NS{z l )-b 1 ^N 1 ^z l i) j>^-.^ V (NS{z 2 )-b^N^z 2 i)G{z l ,z 2 ). (55) 

The reader might be concerned by the slightly unorthodox term in the exponent proportional to N 1 / 3 , however it 
is subdominant to the N S(z) term and its presence does not affect any of the usual arguments of the saddle point 
method; the asymptotic behavior of l|55|) is determined by S(z). From 1(22(1. ((23(1 and (|25f) we see for both the GUE 
and LUE that S(z) has one double saddle point, located at z = — 1. We can deform the contour of integration to a 
new contour which passes through z = — 1 along paths of steepest descent. Note that although in 1(12(1 the integrals 
must be positively oriented, since we have two integrals in 1)55(1 we are free to orient the integrals in the negative 
direction since the consequent minus signs cancel. Let us denote by C the contour consisting of the union of two rays 
of unit length, the first starting at z = e~ lTr / 3 — 1 and ending at z = — 1 and the second starting at z = — 1 and ending 
at z = e"/ 3 - 1. ff we then denote by A any suitable arc such that C U A is a simple closed curve enclosing the origin 
we can write 
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■MO 



.4 



2ttI 



exp 



.4 



dz 2 
2ttI 



exp 



(jVS^-fc 1 / 3 ^ 1 / 3 z 2 e) G( Zl ,z 2 ), 



(56) 



2™ 



T exp 



dz 2 
2ttI 



exp 



(N S(z 2 ) - b l l 3 N 1 ' 3 z 2 G( Zl ,z 2 ) + o(e 2N s ^N~p)- 



(57) 



The error bound in l|57|) holds for all p £E N, and so in what follows we consider p as arbitrarily large. The equality 
between (|57|l and l|56|l can be obtained by noting that we can choose A to consist of two rays lying along the path 
of steepest descent away from the endpoints of C, which extend as far as we like into the right half plane, together 
with an arc to close the contour which we can choose to be as far into the right half plane as desired. With such a 
choice for A one can obtain the required bounds by a straightforward generalization of the usual argument used in 
the saddle point method. For a careful discussion of the saddle point method suitable for this purpose see for example 
Section 2.5 of Ref>i&. We note that Ref>i& refers to the saddle point method as Perron's method. 

Now let us change variables in H57fl according to t = z + 1, so that the vertex of our contour is now at the origin. 
We shall denote the image of C under this change of variables by B. Further, since S(z) is analytic on C we have 



where 



S(t-l)=S(-l) + bj + bj<p{t), 



(58) 



<p(t) 



" gW(-l) 31 tfc _ 3 



fc=4 



S( 3 ) (-1) M 



(59) 



and we have now chosen 



£( 3 >(-l) 
2 ' 
jl, GUE, 
I 2, LUE. 



(60) 
(61) 



We also note that by setting x = 1 in (|31Jl and il55l) we have 




GUE, 

2 - log(2) - in, LUE. 



(62) 



Finally, defining 



F(X,t) :=G(ti -I,i2-l)e 



(63) 



and 



Ef > := / exp ( bN± - £b 1/3 N 1/3 t, 



dU 
2m 



(64) 



we have 



MO 



3 2_/VRc[S(-l)]+2& 1/3 JV 1/3 J 



E 



(S) W (S) 



E 2 °> F(\,t)\ Xz=bNtV3 + o(N~P) 



(65) 



If we now multiply (|65JI by the prefactors required in l(14() using 1)21(1 we obtain 



(WVW 3 ( b 1 ' 3 £ 

o PN 1 + 



2 iV 2 /3 



90 (0 



gi(0 , fla(0 

ATl/3 + AT2/3 



iV 4/3 Ep } EfV(A,t) 



A i= 6iVt?/3 



+ o(iV- p ) 



(66) 



11 




where g m (Q is the coefficient of N m / 3 in the large N fixed £ expansion of 

(67) 

, LUE. 

In 1)66(1 we have presented only terms 0(1/N) in the first factor since this is will be sufficient for our purposes in what 
follows. Higher order terms are easily retained if desired. 

Our work is now essentially done. One expands F(X, t) around t = for fixed A as in l|47|l and then sets Ai = bNt 3 /3, 
analogous to the bulk case. Again, after integration, each value of j in the Maclaurin expansion l|47(l contributes to 
the same order in N. To see this explicitly we can use the following lemma. 

Lemma 1. Let B be the contour consisting of the union of a ray starting at e~ l7r / 3 and ending at the origin, and a 
ray starting at the origin and ending at e n l 3 . For any b > and < (3 < 1/3 we have for large N that 

z m exp ( bN^-^N^z) ^- = (-l)™&-( m+1 V 3 N -(™+i)/3 [M (m)^ + ( e -/»*)] 
\ 3 / 2iri 

where Ai^ m '(£) is the m th derivative of the Airy function Ai(£). 

Proof. This follows from the standard entire contour integral expression for Ai(£) (see e.g^i) by simply changing 
variables z i— » b^^N^^z, and noting that rays defining the contour B can be extended to infinity at the cost of 
introducing exponentially subdominant corrections. □ 

An immediate consequence of Lemma ^ is that with \ = bNtf/3 we have 

( _-i\h+l2+m 1 +m 2 f i \ mi+m.2 

^2/3 E (S) E (S) tT tT X h x i, = J 6(mi+m2+2)/3 Ai^ 1+3 '') (0 Ai^ +3fa ) (0 f^J +0(e-? m ). (68) 

Hence, if we construct the Maclaurin expansion of F(X, t) with A fixed as in l(47|) . and set Ai = bNt 3 /3 and integrate 
using H68fl . we obtain an expansion for 

N^E[ s) E { 2 s) F(X,t)\ Xi=bNt s /3 (69) 
in powers of iV -1 / 3 . Denoting the coefficient of N~ m / 3 in this expansion by c m (£) we have explicitly that 

m + 2 i r Qk Qm+2—k 



^k\{m + 2-k)\ 1212 [dslds 



m+2-k 



F(X,S 1 ,S 2 ] 



Sl ,S 2 =0 



(70) 

\ %= bNtl/3 



The reader might be concerned that according to (|68|l the k\ + k 2 — and fej + fc 2 = 1 terms grow with N; however 
it is not hard to show from l|7l)|l that the coefficients c^ -2 ^^) and c^ -1 -*^) vanish identically. 
We can now express (|66|) in terms of the coefficients c m (£) as 

(bN) 1 / 3 ( b 1 / 3 f \ 1 



o M2/3 ; -wvsy-uvs/ i wivs/-uvs/ i auvsy -xvsyj -M/g 

" ' i /i^ <71) 

+ [•»(£) «(£) + 91 (0 «({) + »({)«(«] + O I - 

We have explicitly displayed terms o(l/N) here, but it straightforward to retain as many terms as desired. The 
expansion constructed from terms o{\/N) however is extremely accurate, as we demonstrate in Figures |3] and ^ and 
it appears to be the numerically optimal order at which to truncate the expansions. 

The explicit forms for the coefficients c m (£) can be constructed from (|70|l and substituted into J7TJ. We can also 
further simplify the Airy derivatives appearing in <|68l) using the Airy differential equation Ai (£) = £ Ai(£) so that 
only Ai(£) and its first derivative appear. We finally obtain the following. 

Proposition 2. Let pn{x) be as defined in ^j). Then with £ fixed, as N — > oo we have the following: 
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FIG. 3: Comparison of the asymptotic expansion %!2\ . shown as the dashed line, and the exact result il l (II . shown as the solid 
line, for the eigenvalue density near the soft edge at x — 1, for the GUE with N = 10. 



For the GUE 



while for the LUE 
(2iV) 1 /3 



PN 1 + 



2iV 2 /3 



(2iV) 2 / 3 



[Ai'(0] 2 -e[Ai(0] 2 
1 



20 
O 



(3£ 2 [Ai(£)] 2 - 2£[Ai'(£)] 2 - 3Ai(0Ai'(0) 



AT2/3 



(72) 



[Ai'(0] 2 -C[Ai(0] 2 + WI [Ai(0] 2 ^ T 



ATi/3 



+ ^(3e 2 [Ai(0] 2 - 2£[Ai'(£)] 2 + (2 - 5a 2 )Ai(0Ai'(C)) ' 



10 

o 



AT2/3 



(73) 



N 



Figures and 01 provide a numerical comparison of the asymptotic expansions given in Proposition[2]with the exact 
results computed using the expression ljl0|l in terms of orthogonal polynomials. 

Note that it appears that the GUE converges much faster than the LUE, since while the two curves in @ are 
almost indistinguishable at N = 10 for the GUE, the asymptotic expansion for the LUE begins to diverge from the 
exact result in already by £ ~ 2, and both curves are rather different from the limiting Airy kernel expression JHJ. 
We investigated the affects of retaining more terms in the expansion for the LUE case; keeping terms 0(1/N) did 
not noticeably change the plots, while keeping terms higher than l/N caused significant divergence of the asymptotic 
expansion from the exact result. The explanation for this is most likely that, as made precise in the next section, 
the edge expansions match onto the bulk expansions, and these in turn become more divergent near the edge at each 
order in l/N. 
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-6 -4 -2 2 4 6 

% 

FIG. 4: Comparison of the asymptotic expansion 1[TA} . shown as the dashed line, and the exact result 11101 . shown as the solid 
line, for the eigenvalue density near the soft edge at x — 1, for the LUE with a = 1/2 and TV = 20. Also shown is the limit as 
N — > oo given by the Airy kernel (JHJ, shown as the dotted line lying below the other two curves. 



V. MATCHING OF THE BULK AND EDGE EXPANSIONS 



In Figures plots of the bulk and edge asymptotic expansions have separately been compared against the exact 
density for N — 10. Although the scale of the independent variable is different, we can see from Figures [3] and 0] that 
the edge asymptotic expansions are accurate approximations to the exact density up to the neighborhood of the first 
local maximum (relative to the edge £ = 0) at least and thus should be used instead of the bulk asymptotic expansion 
in this region. 

At a quantitative level, it is possible to exhibit a matching between the various asymptotic expansions. Suppose in 
(|STT) we set x = 1 + £/2N 2 / 3 , and in (J52JI we set x = 1 + £,/(2N) 2 / 3 , and take £ < and fixed. Expanding the right 
hand sides as an asymptotic series in N gives 



' Vjg cos(4|£| 3 / 2 /3) ' 
7T 2tt|£| 



|g| 3 / 2 cos(4|g| 3 / 2 /3) |£| 3 / 2 sin(4|£| 3 / 2 /3) x 



AT4/3 




1 



20tt J N 2 / 3 

(74) 



q(l + sin(4|C| 3 / 2 /3)) 1 / 1 \ 

WW\ (2Ny/ 3+ °{NV 3 ) (75) 

where the symbol ~ denotes that the asymptotic series have been expanded as specified. An important feature is 
that this procedure mixes the terms which are at different orders in N in l|51J) and (|52(l . 

Let us now compute the £ — > — oo asymptotic expansions of the right hand sides of the first two terms in each of 
(|72|l and J73J), multiplied by 2. Using the fact that for x — > oo (see e.gii£) 

M(-x) ~ _l_cosf^-^ 3 / 2 )-^Lcosf- + -x 3 / 2 )+of4 J N ) 
v ; ./^x 1 ^ V4 3 / 24a; 1 /2 V4 3 / \x 3 / A J 
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we obtain expansions which reproduce the N- independent terms in (|74() and l|75[) . giving furthermore, terms of higher 
order in 1/|£|. In (|75|) the term proportional to l/(2iV) 1 ' 3 is reproduced, and this too is accompanied by terms of 
higher order in 1/|£|. In (|74Jl the term — |£| 3//2 /47riV 2//3 is reproduced, while the other terms proportional to 1/7V 2 / 3 
are out by rational factors. The explanation for the missing higher order terms in 1/|£|, and incorrect rational factors 
is most likely due to the fact that terms of all orders in 1/N in (|51|l and l|52|l contribute to each distinct order in the 
expansions (|74J) and l|75|) . Specifically, from the results exhibited above, it would seem that expanding the complete 
large N asymptotic series for p^ UE (a;) and p^J JE (x) as in 174|) and (|75|l would give precisely the large £ — ► — oo 
expansion of (f72*|l and J73J), extended to all orders in N. 
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